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We have made a calculation of the propagation of positrons from dark-matter particle annihilation in the Galactic 
halo in different models of the dark matter halo distribution using our 3D code, and present fits to our numerical 
propagation Green's functions. We show that the Green's functions are not very sensitive to the dark matter dis- 
tribution for the same local dark matter energy density. We compare our predictions with computed cosmic ray 
positron spectra ("background") for the "conventional" CR nucleon spectrum which matches the local measure- 
ments, and a modified spectrum which respects the limits imposed by measurements of diffuse Galactic 7-rays, 
antiprotons, and positrons. We conclude that significant detection of a dark matter signal requires favourable 
conditions and precise measurements unless the dark matter is clumpy which would produce a stronger signal. 
Although our conclusion qualitatively agrees with that of previous authors, it is based on a more realistic model 
of particle propagation and thus reduces the scope for future speculations. Reliable background evaluation re- 
quires new accurate positron measurements and further developments in modelling production and propagation 
of cosmic ray species in the Galaxy. 

95.35.+d, 98.35.Gi, 98.38.Am, 98.70.Sa 
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I. INTRODUCTION 

Investigations of galaxy rotation, big-bang nucleosynthe- 
sis, and large-scale structure formation imply that a signif- 
icant amount of the mass of the universe consists of non- 
luminous dark matter [fit]. Among the favored particle dark 
matter candidates are so-called weakly interacting massive 
particles (WIMPs), whose existence follows from supersym- 
metric models (see Ref. [|| for a review). If stable, such parti- 
cles could have a significant cosmological abundance at the 
present time. A pair of stable WIMPs can annihilate into 
known particles and antiparticles and it may be possible to 
detect WIMPs in the Galactic halo by the products of their 
annihilations. The difficulty, however, consists in discriminat- 
ing between the products of WIMP annihilation and "back- 
ground" cosmic ray (CR) particles. The smallest background 
arises when considering antiprotons and positrons, secondary 
products of interactions of CR particles with interstellar mat- 
ter, and thus these provide the best opportunity to search for 
dark matter signatures A search for a distinct signature 

in 7-rays from the Galactic halo has also been proposed |Q] . 

Though the microphysics is quite well understood and 
many groups make sophisticated calculations of the spectra of 
annihilation products for numerous WIMP candidates which 
include many decay chains there are still uncertainties 
in the macrophysics which could change the estimated fluxes 
of WIMP annihilation products by 1-2 orders of magnitude, 
making predictions for their detection difficult. In the case 
of antiprotons the strongest evidence would be detection of 
low energy particles (below ~ 1 GeV in the interstellar 
space) but the solar wind and magnetic field sweep low en- 
ergy particles away from the heliosphere, the effect known as 
solar modulation. In the case of 7-rays, a weak signal would 
compete with the flux of Galactic halo 7-rays, the uncertain 
background of extragalactic photons, and a contribution from 
unresolved sources. 



The most promising is perhaps the positron signal since it 
can appear at high energies where the solar modulation is neg- 
ligible, but its strength depends on many details of propaga- 
tion in the Galaxy. The "leaky box" model is often used [pra], 
a simplified approach which may not be applicable in the case 
of positrons. The most accurate propagation model applied 
so far, the diffusion model is analytical and thus is sub- 
ject to certain simplifications, e.g. the positron source func- 
tion is treated as being dependent only on the radial cylindri- 
cal coordinate R, the assumption of spatially uniform inter- 
stellar radiation and magnetic fields, and some other minor 
details. On the other hand, progress in CR positron mea- 
surements is anticipated since several missions operating or 
under construction are capable of measuring positron fluxes 
up to 100 GeV (e.g. gas-RICH/CAPRICE and PAMELA ex- 
periments [^]). Therefore, more accurate calculation of the 
positron propagation is desirable. 

We have developed a numerical method and corresponding 
computer code (GALPROP)[] for the calculation of Galactic 
CR propagation in 3D [jiot|. The rationale for our approach 
was given previously Jl0|-|15]]. Briefly, the idea is to develop 
a model which simultaneously reproduces observational data 
of many kinds related to cosmic -ray origin and propagation: 
directly via measurements of nuclei, antiprotons, electrons, 
and positrons, indirectly via 7-rays and synchrotron radia- 
tion. These data provide many independent constraints on any 
model and our approach is able to take advantage of this since 
it aims to be consistent with many types of observation. 

The code is sufficiently general that new physical effects 
can be introduced as required. Its capability includes primary 
and secondary nucleons, primary and secondary electrons, 



*Our model including software and datasets is ava ilable at 
nttp://www.gamma.mpe-garching.mpg.de/^aws/aws.html 



secondary positrons and antiprotons, as well as 7-rays and 
synchrotron radiation. The basic spatial propagation mecha- 
nisms are diffusion and convection, while in momentum space 
energy loss and diffusive reacceleration are treated. Fragmen- 
tation, secondary particle production, and energy losses are 
computed using realistic distributions for the interstellar gas 
and radiation fields. We aim for a "standard model" which 
can be improved with new astrophysical input and additional 
observational constraints. 

In this paper we use our model for calculation of positron 
propagation in different models of the dark matter halo distri- 
bution. We compare our predictions with evaluated cosmic 
ray positron spectra ("background") for the "conventional" 
CR nucleon spectrum which matches the local measurements, 
and for a modified spectrum which respects the limits imposed 
by measurements of diffuse Galactic 7-rays, antiprotons, and 
positrons. To be specific we will further discuss neutralino 
dark matter, although our results can be easily adopted for any 
other particle dark matter candidate. 



II. BASIC FEATURES OF THE GALPROP MODELS 

The GALPROP models have been described in full detail 
elsewhere [|l(J; here we just summarize briefly their basic fea- 
tures. 

The models are three dimensional with cylindrical symme- 
try in the Galaxy, and the basic coordinates are (R, z,p) where 
R is Galactocentric radius, z is the distance from the Galactic 
plane and p is the total particle momentum. In the models the 
propagation region is bounded by R — Rh, z = ±z^ beyond 
which free escape is assumed. 

The propagation equation we use for all CR species is writ- 
ten in the form: 

01b —* —* —* 3 3 I 
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where ip = ip(f,p,t) is the density per unit of total parti- 
cle momentum, ip(p)dp = 4irp 2 f(p) in terms of phase-space 
density f(p), q{r,p) is the source term, D xx is the spatial 
diffusion coefficient, V is the convection velocity, reacceler- 
ation is described as diffusion in momentum space and is de- 
termined by the coefficient D pp , p = dp/dtis the momentum 
loss rate, r/ is the time scale for fragmentation, and r r is the 
time scale for the radioactive decay. 

An assumption is free escape of particles at the halo bound- 
aries. Under certain simplifications it translates into the re- 
quirement that the number density of particles at the bound- 



aries is zero: 



ip(Rh,z,p) = ip(R, ±z h ,p) = 0. 



(2) 



This is an approximation, but since the number density of 
particles there is presumably small, it is resonable to assume 



that this should not affect much the particle distribution in the 
Galaxy. 

The numerical solution of the transport equation 
is based on a Crank-Nicholson JT^ ] implicit second-order 
scheme. Since we have a 3-dimensional (R, z,p) problem we 
use "operator splitting" to handle the implicit solution. We ap- 
ply the implicit updating scheme alternately for the operator 
in each dimension in turn, keeping the other two coordinates 
fixed. A check for convergence is performed by computing 
the timescale ijj / (dijj / dt) from eq. (lj) and requiring that this 
be large compared to all diffusive and energy loss timescales. 
The details of our method are fully explaned in JTo| ] . 

For a given zh the diffusion coefficient as a function of mo- 
mentum and the reacceleration parameters is determined by 
CR Boron-to-Carbon (B/C) ratio. Reacceleration provides a 
natural mechanism to reproduce the B/C ratio without an ad- 
hoc form for the diffusion coefficient. The spatial diffusion 
coefficient is taken as f3Do(p/ po) s . Our reacceleration treat- 
ment assumes a Kolmogorov spectrum with 5 = 1/3. For the 
case of reacceleration the momentum-space diffusion coeffi- 



cient Dpp is related to the spatial coefficient D xx 1170. The 
injection spectrum of nucleons is assumed to be a power law 
in momentum, dq(p) /dp oc p~ 7 for the injected particle den- 
sity, if necessary with a break. 

The total magnetic field is assumed to have the form 



Btot = B e 



-(R-R @ )/R B -\z\/zi 



(3) 



The values of the parameters (Bo, Rb, Zb) are adjusted to 
match the 408 MHz synchrotron longitude and latitude dis- 
tributions. The interstellar hydrogen distribution uses HI and 
CO surveys and information on the ionized component; the 
Helium fraction of the gas is taken as 0.11 by number. En- 
ergy losses for electrons by ionization, Coulomb interactions, 
bremsstrahlung, inverse Compton, and synchrotron are in- 
cluded, and for nucleons by ionization and Coulomb inter- 
actions. The distribution of cosmic-ray sources is chosen to 
reproduce the cosmic-ray distribution determined by analysis 
of EGRET 7-ray data Q and was described in Ref. 0. 

Positron production is computed as described in Ref. [[TT|]; 
this includes a critical reevaluation of the secondary decay cal- 
culations. 

Gas related 7-ray intensities are computed from the emis- 
sivities as a function of (R, z, E 7 ) using the column densities 
of HI and H2 for Galactocentric annuli based on 21 -cm and 
CO surveys jl8| ] . The interstellar radiation field (ISRF), which 
is used for calculation of the inverse Compton (IC) emission 
and electron energy losses, is calculated based on stellar pop- 
ulation models and COBE results, plus the cosmic microwave 
background. Our results for diffuse continuum 7-rays, syn- 
chrotron radiation, and a new evaluation of the ISRF are given 
in Ref. [£2|]. 

An overview of our results is presented in Ref. [ pj| ] and 
full results for protons, Helium, positrons, and electrons in 
Ref. 0. The evaluation of the B/C and 10 Be/ 9 Be ratios, 
evaluation of diffusion/convection and reacceleration models, 
and full details of the numerical method are given in Ref. JTo|]. 
Antiprotons have been evaluated in the context of the "hard 
interstellar nucleon spectrum" hypothesis in Ref. [p4[]. 
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TABLE I. Parameters of the dark matter profiles [N9h. 



Model 


po, GeV cm 3 


T c , kpc 


"isothermal" 


0.43 


2.8 


Evans 


0.51 


7.0 


alternative 


0.38 


0.9 



III. GREEN'S FUNCTIONS 

To make a prediction of the positron flux at the solar po- 
sition, one needs to know the source function /(e) which de- 
scribes the spectrum of positrons from neutralino annihilation, 
and the Green's function G(E, e) for their propagation in the 
Galaxy. Then the positron flux is a convolution 



dF 
dE 



(av. 



pl 



J deG(E,e)Y,B l f l (e) 



[cm~ 2 s _1 sr -1 GeV" 1 ], (4) 



where (av) is the thermally averaged annihilation cross sec- 
tion, po is the local dark matter mass density, m x is the neu- 
tralino mass, Bi is the branching ratio into a given final state 
i. The Green's function thus includes all details of the dark 
matter mass distribution and Galactic structure (diffusion co- 
efficient, spatially and energy dependent energy losses etc.). 

Following Kamionkowski and Kinkhabwala [|l9[| we con- 
sider three different dark matter mass density profiles which 
match the Galactic rotation curve. The canonical "isothermal" 
sphere profile, 



R% 



(5) 



where r c is the core radius, Rq = 8.5 kpc is the solar distance 
from the Galactic center, and r 2 = R 2 + z 2 is the spherical 
radial coordinate. The spherical Evans model, 



P(r) = Po 



(r 2 +R 2 ) 2 irl+r 2 
3r 2 + R% (r 2 +r 2 ) 2 ' 



and an alternative form which also might be empirically ac- 
ceptable, 



p(r) = po 



(r c + r) 2 



(7) 



Note that po and r c for each model must be fitted to the ro- 
tation curve, and therefore they are different for each model 
(see Table |). These profiles are plotted in Fig. [j]. 

For each given model we calculate the Green's function 
G(E, e) defined in Eq. (Q), which gives the positron flux at 
the solar position corresponding to the positron source func- 
tion in the form of a Dirac (5-function in energy. The positron 
propagation is calculated in a model which was tuned to match 
many available astrophysical data [ |io[|l2| ]. Since the halo size 
in the range zu = 4 — 10 kpc is favored by our analyses of 
B/C and 10 Be/ 9 Be ratios and diffuse Galactic 7-ray emission 
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FIG. 1. The radial profiles of the spherical halo models: the 
canonical isothermal model (solid line), Evans model (long dashes), 
alternative model (dots), and uniform distribution (p = 0.4 GeV 
cm -3 , short dashes). 



[ }10Ul2p , we consider two cases z% = 4 and 10 kpc which 
provide us with an idea of the possible limits. The preferred 
neutralino mass range following from accelerator and astro- 
physical constraints is 50 GeV < m x < 600 GeV [p0|], and 
we consider positron energies e < 824 GeV which cover this 
range. 

For the case of a uniform dark matter mass distribution 
P( r ) = (p) = const we compare our results with simple 
analytical Green's functions for the leaky-box model for two 
cases, where the positron containment time is a constant pa- 
rameter, t = T , 



G t {E,e) 



1 



4tt£ E 2 



exp 



T0t 



(e - E), (8) 



(g) and when it varies with energy, r(e) = rj/e, 



G 2 (E,e) 



c 1 
I^e 2 



(9) 



where c is the speed of light, n is a constant, £ is the energy 
loss constant de/dt = £e 2 , and 6(x) is the Heaviside step 
function. 

Fig. |2] shows functions E 2 G(E, e) calculated in our model 
for Galactic halo sizes zn = 4 and 10 kpc compared to 
the leaky-box functions Gx^{E, e). For this comparison we 
adopted the following parameters [Q]: G\. £ = 1.11 x 10 
yr- 1 GeV~\ r = 10 7 yr; G 2 : £ = 1.52 x 10" 9 yr" 1 GeV 
i] = 2 x 10 8 yr GeV. It is clear that the leaky-box model does 
not work here, moreover a resonable fit to our G-functions is 
impossible for any combination of £ and To (or 77). The dif- 
ference in the normalization at maximum (E = e) is mainly 
connected with our accurate calculation of the ISRF which is 
responsible for the energy losses. 
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FIG. 2. Calculated G-functions for the uniform dark matter dis- 
tribution, z h = 4 kpc and 10 kpc, for e = 25.76, 103.0, 412.1 
GeV (solid lines). The leaky-box functions Gi and G2 are shown 
by dashed and dotted lines respectively. The units of the abscissa are 
10 25 GeV cm sr _1 . 



Fig. |3| shows our calculated G-functions for different mod- 
els of the dark matter distribution: "isothermal", Evans, and 
alternative. The curves are shown for two halo sizes Zh = 4 
and 10 kpc and several energies e = 1.03, 2.06, 5.15, 10.3, 
25.8, 51.5, 103.0, 206.1, 412.1, 824.3 GeV. At high energies, 
increasing positron energy losses due to the IC scattering com- 
pete with the increasing diffusion coefficient, while at low en- 
ergies increasingenergy losses due to the Coulomb scattering 
and ionization JlC| ] compete with energy gain due to reaccel- 
eration. The first effect leads to a smaller sensivity to the halo 
size at high energies. The second one becomes visible below 
~ 5 GeV and is responsible for the appearance of accelerated 
particles with E > e. 

It is interesting to note that for a given initial positron en- 
ergy all three dark matter distributions provide very similar 
values for the maximum of the G-function (on the E 2 G(E, e) 
scale), while their low-energy tails are different. This is 
a natural consequence of the large positron energy losses. 
Positrons contributing to the maximum of the G-function 
originate in the solar neighbourhood, where all models give 
the same dark matter mass density [see eq. (Q) for the def- 
inition of the G-function]. The central mass density in these 
models is very different (Fig. [j]), and therefore the shape of the 
tail is also different since it is produced by positrons originat- 
ing in distant regions. As compared to the isothermal model, 
the Evans model produces sharper tails, while the alternative 
model gives more positrons in the low-energy tail. At inter- 
mediate energies (~ 10 GeV) where the energy losses are 
minimal, the difference between z;, = 4 and 10 kpc is max- 
imal. Also at these energies positrons from dark matter par- 
ticle annihilations in the Galactic center can contribute to the 
predicted flux. This is clearly seen in the case of the alterna- 
tive model with its very large central mass density (Fig. ||c, 
z h = 10 kpc). 

To provide the Green's function for an arbitrary positron 
energy, which is necessary for prediction of positron fluxes in 
the case of continuum positron source functions (as will be 



required if one considers secondary, tertiary etc. decay prod- 
ucts), we made a fit to our numerical results. Since a reson- 
able fit using the leaky-box Green's functions is impossible 
we have chosen the function 



G(E,e) 



10 2 



E 2 



| 10 alog 



(e-E) 



^Q^iog^+ziogB+y^^j [ cmsr -i GeV" 1 ], (10) 

which allows us to fit our numerical functions with accu- 
racy better than 10% over a decade in magnitude (on the 
E 2 G(E, e) scale). Here the first term fits the low energy 
tail, the second term fits the right-hand-side part of the G- 
functions and represents the effect of reacceleration, E is in 
GeV, and a(e), 6(e), c(e), w(e), x(e), y(e) are the fitting pa- 
rameters. Though a better fit is possible by using more com- 
plicated functions, we try to minimize the number of fitting 
parameters while still providing resonable accuracy. Besides, 
the accuracy of our propagation model is not better than 10%, 
being limited by the accuracy of the astrophysical data input. 
The numerical values of the fitting parameters are given in 
Tables ^1 - III for the three models discussed. At intermediate 
energies the parameters can be interpolated. The cubic spline 
(or square spline for w, x, y) provides ~ 10% accuracy for 
the G-functions when interpolating on the logarithmic energy 
scale. 



IV. POSITRON FLUXES 
A. Positrons from the dark-matter particles annihilation 

When neutralinos annihilate in the Galactic halo they pro- 
duce quarks, gluons, leptons, and other particles which via 
hadronization and/or decays give rise to secondary positrons. 
One can expect to get both monoenergetic positrons (en- 
ergy m x ) from direct annihilation into e + e~ and continuum 
positrons from the other annihilation channels. In general, the 
direct e + e~ annihilation channel is severely suppressed with 
a branching ratio of order 10~ 5 [Bh, though some classes of 
models allow a larger branching ratio to be obtained. Also in 
some cases, e.g. if annihilation occurs near a pole in the cross 
section, (av) can be quite large [ 2 1 1 which can compensate for 
a small branching ratio. The (quasi-) monoenergetic positron 
line, if it is strong enough, is easier to identify. In contrast, the 
hadronization and/or decay cascades lead to the appearance of 
"continuum" positrons with correspondingly degraded energy 
thus making worse the signal/background ratio. 

For a neutralino heavier then W (Z°) boson, annihilation 
to W ± - or Z°-boson pairs will be significant followed by the 
direct decay of M /+ 's and Z°'s to lepton pairs where the direct 
positron channel accounts for 11% and 3.4% of W + and Z Q 
decays, respectively. In the case where the neutralino is a pure 
Higgsino state, the annihilation cross section to W - or Z°- 
boson pair increases rapidly above the threshold and reaches 
a maximum of (av) w 3 x 10~ 25 cm 3 s _1 (for W + W~) 
and w 2 x 10 



25 cm 3 s- 1 (for Z°Z°) at about 110 GeV and 
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E A 2*G 'Isothermal' dark matter distribution 




E, GeV 

FIG. 3. Calculated G-functions for different models of the dark matter distribution: (a) "isothermal", (b) Evans, (c) alternative. Upper 
curves z h = 10 kpc, lower curves z h = 4 kpc, e = 1.03, 2.06, 5.15, 10.3, 25.8, 51.5, 103.0, 206.1, 412.1, 824.3 GeV. The units of the 
abscissa are 10 25 GeV cm sr 



5 



TABLE II. Fitting parameters of the Green's functions. 



"isothermal" model Evans model alternative model 



e, GeV 


a 


b 


c 


a 


b 


c 


a 


b 


c 


z h = 4 kpc 




















1.03 


-1.9732 


2.3448 


-0.1340 


-1.9937 


2.6865 


-0.2421 


-2.0146 


1.8984 


0.0364 


2.06 


-2.4853 


3.2517 


-0.6564 


-2.2043 


3.5919 


-0.8762 


-2.9275 


2.7980 


-0.3409 


5.15 


-2.6365 


4.4743 


-1.6189 


-2.4577 


4.7878 


-1.9748 


-3.0137 


4.1174 


-1.1397 


10.30 


-1.9555 


4.4101 


-2.2099 


-1.9686 


4.8817 


-2.6923 


-2.0712 


3.9274 


-1.6167 


25.76 


-1.1684 


3.7535 


-2.6853 


-1.3763 


4.5614 


-3.4205 


-0.9412 


2.7981 


-1.8075 


51.52 


-0.8469 


3.3985 


-3.0180 


-1.1035 


4.3731 


-3.9446 


-0.4843 


2.0882 


-1.8434 


103.00 


-0.6979 


3.3043 


-3.4776 


-0.9532 


4.3484 


-4.5560 


-0.2718 


1.6991 


-1.9641 


206.10 


-0.6173 


3.3187 


-4.0118 


-0.8078 


4.1933 


-5.0264 


-0.2586 


1.8175 


-2.4451 


412.10 


-0.5337 


3.2338 


-4.4403 


-0.6719 


3.9382 


-5.3478 


-0.3292 


2.2664 


-3.2998 


824.30 


-0.4255 


2.9458 


-4.6095 


-0.5078 


3.4153 


-5.2860 


-0.3330 


2.4540 


-3.9563 


Zh = 10 kpc 




















1.03 


-2.9531 


2.2688 


0.2037 


-2.2651 


2.7232 


0.0540 


-3.9088 


1.6600 


0.4355 


2.06 


-4.3610 


4.1399 


-0.4085 


-3.3310 


4.1276 


-0.6391 


-5.3710 


4.0111 


-0.0642 


5.15 


-4.6245 


6.7914 


-1.9188 


-3.7297 


6.2467 


-2.0800 


-5.8784 


7.5056 


-1.6622 


10.30 


-3.6318 


7.1054 


-2.9096 


-3.0093 


6.5384 


-3.0367 


-4.7439 


8.2050 


-2.8244 


25.76 


-2.3450 


6.2626 


-3.6392 


-2.0541 


5.9838 


-3.8567 


-2.8799 


6.8568 


-3.4165 


51.52 


-1.6556 


5.3936 


-3.8732 


-1.5142 


5.3295 


-4.1955 


-1.8345 


5.3634 


-3.3196 


103.00 


-1.1707 


4.5572 


-3.9405 


-1.1606 


4.8040 


-4.4885 


-1.1501 


4.0735 


-3.0763 


206.10 


-0.8276 


3.8034 


-3.9013 


-0.8889 


4.2660 


-4.6500 


-0.7058 


3.0354 


-2.7952 


412.10 


-0.5822 


3.1489 


-3.8069 


-0.6680 


3.6993 


-4.6643 


-0.4132 


2.2082 


-2.5103 


824.30 


-0.3984 


2.5623 


-3.6580 


-0.4917 


3.1504 


-4.5846 


-0.2256 


1.5851 


-2.2758 



TABLE III. Fitting parameters of the Green's functions. 





"isothermal" model 






Evans model 




alternative model 




e, GeV 


w 


X 


y 


w 


X 


y 


w 


X 


V 


z h = 4 kpc 




















1.03 


-4.9292 


-0.8786 


-0.1115 


-4.0224 


-\Mn 


-0.2053 


-5.9792 


-0.1569 


0.0412 


2.06 


-7.2475 


1.0942 


0.4742 


-6.5523 


0.1919 


0.6011 


-8.1532 


2.2559 


0.3291 


5.15 


-8.9618 


2.1785 


3.1988 


-8.0219 


0.4044 


3.9380 


-9.3685 


3.1591 


2.7479 


Zh = 10 kpc 




















1.03 


-4.3201 


-1.1227 


0.2324 


-3.5986 


-1.4890 


0.0867 


-5.3468 


-0.6192 


0.4554 


2.06 


-6.0920 


0.2541 


0.9683 


-4.2529 


-1.4473 


1.1879 


-8.0942 


2.1474 


0.7739 


5.15 


-6.5457 


-1.1929 


4.7067 


-10.0800 


3.5947 


3.0000 


-9.0223 


3.1302 


3.0000 



120 GeV, respectively [Q]. For unpolarized W + -boson the 
decay is isotropic in the rest frame, which produces a uniform 
positron distribution in the laboratory system: 



m x Pw 



(11) 



where f3w is the iy + -boson speed in the laboratory system, 
ande± = ^m x (l±Pw)- Since the CR positron spectrum falls 
as ~ E~ 3 3 above several GeV [|22"|], and the signal strength is 
proportional to e~ 2 m~ 3 , the signal/background ratio is maxi- 



mal near m-t 



B. Positron "background" 

An important issue in interpretation of the positron mea- 
surements is evaluation of the "background", positrons arising 



from CR particle interactions with interstellar matter. Though 
the parameters of the propagation and the Galactic halo size 
can be fixed in a self-consistent way using CR isotope ratios, 
the ambient CR proton spectrum on the Galactic scale remains 
quite uncertain. 

The only possibility to trace the spectrum of nucleons on 
a large scale is to observe secondary products such as diffuse 
7-rays, positrons, and antiprotons. The EGRET data show en- 
hanced 7-ray emission above 1 GeV in comparison with cal- 
culations based on locally measured ("conventional") proton 
and electron spectra [p3j]. This can be interpreted as imply- 
ing that the average spectra of particles in the Galaxy can dif- 
fer from what we measure locally, due to details of Galactic 
structure and, in the case of electrons, large energy losses. A 
possible solution could be a hard interstellar proton spectrum 
p4], or an electron spectrum which is on average harder than 
that locally observed [ |25| , [26[ ] due to the spatially inhomoge- 
neous source distribution and energy losses. 
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The first possibility has been studied in detail in relation to 
antiprotons and positrons [ |l2|Jl4| ], showing that the resulting 
particle fluxes are too large. Taken together, the antiproton 
and positron data provide rather substantial evidence against 
the idea of explaining the > 1 Ge V 7-ray excess with a hard 
nucleon spectrum. 

The second possibility looks more plausible and detailed 
calculations [ |l2pq ] showed that the 7-ray excess could in- 
deed be explained in terms of IC emission from a hard electron 
spectrum, but the fit to the EGRET spectral shape is still not 
very good. It can be further improved by allowing some free- 
dom in the nucleon spectrum at low energies, which is possi- 
ble since solar modulation affects direct measurements of nu- 
cleons below 20 GeV. Because of the hard electron spectrum 
the required modification to the nucleon spectrum is moder- 
ate; as expected the predictions for antiproton and positron 
fluxes are larger than in the "conventional" model but still 
within the allowed limits [|l2|] . 

In order to show the effect of varying of the ambient pro- 
ton spectrum, we compare our results with two models for the 
CR positron "background". These are a "conventional" model 
(model C) which reproduces the local directly measured pro- 
ton and Helium spectra above 10 GeV (where solar modula- 
tion is small), and a model with modified nucleon spectrum 
(model HEMN), which is flatter below 20 GeV and steeper 
above, and arises from our analysis of Galactic diffuse 7-ray 
emission. The "background" spectra are slightly dependent 
on the halo size. Since all secondary particles are produced 
in the Galactic plane, increasing the halo size results only in 
a small decrease of the flux at high energies due to larger en- 
ergy losses. Th e propa gation parameters for these models are 
given in Refs. [ |lQ , |l2| ], and the formalism for calculation of 
secondary positrons is described in Ref . [pTJ| . 



0.0001 



C. Calculations 

We do not intend to make sophisticated calculations of 
positron spectra resulting from numerous decay chains such 
as best done by, e.g., Baltz and Edsjo [|[] for many WIMP 
candidates. Instead, for illustration purposes, we simplify 
our analysis by treating the annihilation to W and Z°-pairs. 
For m x < mw we consider only the direct annihilation to 
e + e~ pairs. In the first case we use the cross sections for a 
pure Higgsino [||] and the production source function given 
by Eq. ([fj]), in the latter case we take B ■ (av) = 3 x 10~ 28 
cm 3 s _1 and monoenergetic positrons. These parameters can 
be considered as optimistic, but possible [^|^]. To maximize 
the signal we further choose the Galactic halo size as 10 kpc. 

Fig. |] shows our predictions for the two CR positron "back- 
ground" models together with HEAT data [^2] and positrons 
from neutralino annihilation. It is seen that the predicted sig- 
nal/background ratio has a maximum near m x ~ m-w, while 
even in the "conventional" model the background is nearly 
equal to the signal at its maximum. It is however interesting 
to note that our calculations in this model show some excess 
in low energy (< 10 GeV) positrons where the measurements 




1e-08 



10 100 
E, GeV 

FIG. 4. Our predictions for two CR positron "background" mod- 
els (C and HEMN: heavy solid lines), and positron signals from neu- 
tralino annihilation for m x = 5.15, 10.3, 25.8, 103.0, 206.1, 412.1 
(thin solid lines): (a) zh — 4 kpc, (b) Zh = 10 kpc. In the case of 
m x = 103.0 GeV, the signal plus background (model C) is shown 
by the dotted line. Data and the best fit to the data (dashes) are from 
Ref. Eh (HEAT collaboration). 



are rather precise but the solar modulation is also essential. 
If this excess testifies to a corresponding excess in interstel- 
lar space and if the positron background correspond to our 
"conventional" calculations, it could be a hint for the pres- 
ence of dark matter |||7|]. Our HEMN model fits the HEAT 
data better (no excess) and thus provides more background 
positrons. (This shows that in principle a good fit to positron 
data, which is consistent also with other measurements such 
as 7-rays and antiprotons is possible without any additional 
positron source.) Under such circumstances a significant de- 
tection of a weak signal would require favourable conditions 
and precise measurements. Though this our conclusion qual- 
itatively agrees with that of Baltz and Edsjo [|p and several 
earlier papers, it is based on a more realistic model of parti- 
cle propagation and thus reduces the scope for future specula- 
tions. 

We should mention, however, a possibility which could 
increase the signal by orders of magnitude. Relatively 
small fluctuations in the dark matter density distribution will 
strongly increase the positron yield; this would be the case if 
the dark matter halo is clumpy [[28|]. But even if such a signal 
is detected, its correct interpretation will require reliable back- 
ground calculations and thus emphasizes the necessity for fur- 
ther developments in modelling production and propagation 
of CR species in the Galaxy. 



V. CONCLUSION 

We have made a calculation of propagation of positrons 
from dark-matter particle annihilation in the Galactic halo us- 
ing our 3D model which aims to reproduce simultaneously 
observational data of many kinds related to cosmic-ray ori- 
gin and propagation: directly via measurements of nuclei, an- 
tiprotons, electrons, and positrons, indirectly via 7-rays and 
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synchrotron radiation. We use this model for the calculation 
of positron propagation in different models of the dark matter 
halo distribution and present fits to our numerical propagation 
Green's functions. 

We have shown that the Green's functions are not very 
sensitive to the dark matter distribution for the same local 
dark matter energy density. This is a natural consequence of 
the large positron energy losses. The differences in the cen- 
tral dark matter mass density lead to different shapes of the 
Green's function low-energy tail, since this involves positrons 
originating in distant regions. As compared to the isothermal 
model, the Evans model produces sharper tails, while alterna- 
tive model gives more positrons in the low-energy tail. 

We compare our predictions with the computed CR 
positron "background" for two models of the CR nucleon 
spectrum. We conclude that a significant detection of a dark 
matter signal requires favourable conditions and precise mea- 
surements unless the dark matter is clumpy which would pro- 
duce a stronger signal. Though our result qualitatively agrees 
with that of previous authors, it is based on a more realistic 
model of particle propagation and thus provides a firmer basis 
for conclusions. 

A correct interpretation of positron measurements requires 
reliable background calculations and thus emphasizes the ne- 
cessity for further developments in modelling production and 
propagation of CR species in the Galaxy. The ambient proton 
spectrum is of primary importance; its study requires com- 
bined approach which allows to evaluate 7-ray, antiproton and 
other data simultaneously. 
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